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Abstract 



We report a detailed experimental and theoretical investigation of the effect 



^ ■ of residual strain, and strain relaxation, which manifests itself at the Si/Si02 

)£? • interfaces in commercial silicon-on-insulator (SOI) wafers. SOI material is 

^^ ■ made of a single-crystal silicon overlayer (SOL) on top of an insulator (buried 

o: 

^D . Si02 layer) sitting on a handle silicon wafer. Infrared reflectivity spectra show 



that the buried Si02 layer relaxes continuously when thinning the SOL. At the 



i-rt ' same time the SOL surface roughness and the linewidth of optical phonons in 

a: 

O ■ Si near the Si/Si02 interface (probed by micro-Raman specroscopy) increase. 



In the as-delivered wafers, this comes from a slight expansion of Si on both 
sides of the buried Si02 layer which, conversely, is compressed. Thinning 
the SOL modifies these initial equilibrium conditions. To get quantitative 
results, we have modeled all our Raman spectra using a theory of inhomo- 
geneous shift and broadening for optical phonons, which takes into account 
the phonon interaction with the static strain fluctuations. From the variation 
of linewidth versus interface distance, we have found that the mean squared 



strain continues to relax in the bulk of the wafer through a depth on the 

order of several fim. We also show that the SOL surface roughness is related 

to strain fluctuations near the Si/Si02 interfaces. 
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I. INTRODUCTION 

Silicon-On-Insulator (SOI) technology emerged in the early 80s and became increasingly 
popular [|l|, [0. Initially limited to space and military applications where radiation hardness 
is more important than cost and performance, SOI material was rapidly considered for 
PD (Partially Depleted) and FD (Fully Depleted) C-MOS (Complementary-Metal Oxide 
Semiconductor) technology. The net consequence is that IBM announced recently mass- 
production of ICs (Integrated Circuits) on SOI using PD-MOS architectures |^. 

The main reason for success in the case of ICs on SOI is that unlike conventional bulk 
silicon they do not have any direct electrical connection between the active devices and the 
underlying substrate. This renders easier the production of low power / low voltage / high 
frequency ICs. However, since FD-MOS remains difficult to master, there is still room for 
deep scientific interest. In this work, we report a detailed consideration of the strain and 
strain relaxation effects which appear in such non-traditional bonded silicon substrates. 

Whatever the manufacturing technology, all standard SOI wafers have the appearance 
shown in Fig. 1. On top is a thin film of single-crystalline silicon, next a buried oxide film 
and, finally, a thick (handle) silicon wafer. Usually the Silicon-OverLayer (SOL) is 0.2 /im- 
thick and the Buried-OXide (BOX) is 0.4 /xm-thick. Basically two different techniques can 
be used to produce the buried oxide layer. One is oxygen-ion implantation, which results 
in SIMOX (Separation by IMplantation of Oxygen). The second is hydrophilic bonding of 
two (previously oxidized) Si wafers, which results in BESOI (Bond and Etch back SOI) and 
Unibond [Q. Of course, the topmost SOL and the BOX should be defect-free, stress- free, 
and uniform in thickness. The handle Si wafer should be also defect-free. 

The real situation departs from this ideal viewpoint. From a recent survey of literature 
data, the SOL surface roughness ranges from less than 2 A in Unibond to 5A in high dose 
SIMOX, whereas the density of dislocations in the SOL runs from 100 cm~^ in Unibond to 
10^ cm~^ in SIMOX. Since any direct bonding of two oxidized wafers (or any oxygen-ion 
implantation followed by a high temperature annealing step) must introduce a finite amount 



of stress at the interface 0, it is obvious that the primary source of defects in the SOL is the 
stress located at the Si/Si02 interface. Upon relaxation, this stress leads to the generation 
of dislocations which, in turn, affect the silicon parts of the material. 

The study of the interaction between dislocations and a strain field is an old problem. 
In the recent years, it has been actively pursued (see, for instance, Refs. [§[ - |§|), but many 
qualitative observations are still not satisfactorily understood. For instance, in recent work 
1^, we have reported some preliminary results concerning the strain in SOI, which need to be 
quantitatively discussed. We have used SIMOX and Unibond materials and, in both cases, 
started from ~ 200 nm (nominal) SOL thickness. Then we thinned progressively down to 
~ 5 nm and observed three different phenomena: (i) a change in the 1050 cm~^ infrared 
reflectivity structure associated with the internal vibrational modes of the Si04 tetrahedra 
located in the BOX, (ii) an increase in the SOL surface roughness probed by AFM (Atomic 
Force Microscopy) and, (iii) a variation in the linewidths of the optical phonons in the bulk 
of the handle Si-wafer as a function of the distance from the BOX/wafer interface. This 
was probed using micro-Raman spectroscopy. It was found that roughly speaking, given a 
specific SOL thickness, the linewidth becomes narrower as the distance from the interface 
increases. Conversely, given a distance, the linewidth was larger for the samples with the 
thinner overlayer. This behavior corresponds with surface roughness modification. The 
roughness of the free surface increases when the thickness of the overlayer decreases. 

Since the same distance-to-interface dependence of the linewidth is found in different 
kinds of SOI (BESOI and SIMOX materials) it is necessary to propose a common model to 
explain the physical origin of this phenomenon. For instance, because the natural width of 
a Raman phonon line is mainly determined by the anharmonic decay of an optical phonon 
(with frequency ujq) into two phonons of opposite wave vectors, additional contributions near 
the interface could result from the decay into two surface phonons or, in the long-wave limit, 
into two acoustic surface Rayleigh waves. The excited surface phonons have frequencies ujo/2 
and wave vectors k = ujQ/2stC (where Sj is the transverse-acoustic velocity and the constant 
(, of the order of 0.9, depends on the elasticity constants of the material). Therefore, in 



the direction normal to the surface, the corresponding penetration depth S ~ l/Zc, being 
on the order of wavelength, is not larger than a few tens of lattice parameters a. Then the 
contribution to the Raman spectra must be of the order of 6/ds, where dg is the diameter 
of the laser spot. Such a small contribution is not observable in our experiments. 

An alternative reason for broadening could be related to imperfections localized near the 
interface. In a single crystal only the phonon modes with a well defined wave-vector (close 
to the zone-center because the incident and scattered photons have very small momenta) 
can contribute to the Raman spectra. All departures from this "perfect crystal" conditions 
are known as a "relaxation of the momentum selection rule" |jTO|. A simple (qualitative) de- 
scription is usually obtained ||TT| , by integrating the natural (Lorentzian) phonon line shape 
with a Gaussian distribution function of phonon states (the width of which is determined 
by the space scale of the imperfections). However, such a description does not explain the 
physical source of broadening. 

In this work we challenge a model which takes quantitatively into account the scattering 
of phonons by structural imperfections. The central point is strain induced at the Si/Si02 
interface because of the large difference in thermal expansion coefficients. This results in 
structural defects (or misfit dislocations) in the adjacent layers. Away from the interface 
the strain relaxes at some finite distance, involving new threading dislocations. The effect 
of the smooth part of the strain on the phonon frequencies can be compared with the 
infiuence of pressure. With this respect, our theory is consistent with the well known results 
for the pressure dependence of phonon frequencies. This was previously observed under 
homogeneous stress in Si, Ge and GaAS [|12| or thin films of cubic SiC on Si [|1^, |T^ and 



results in well-known phonon shift and splitting. 

The specific effect of strain fluctuations is to give rise to a phonon scattering process 
which induces additional phonon width and shift. Despite the fact that the effect of disorder 
on vibrational modes has been studied for a long time |^, a detailed model has only been 
proposed recently. Taking into account the optical phonon scattering due to disorder, the 
theory describes the inhomogeneous phonon width associated with phonon scattering by the 



strain fluctuations and explains satisfactorily all relevant experimental results ||T^ . 



One important point to notice is that the most distinctive features of the optical-phonon 
scattering by static imperfections is the influence of the DOS singularity at the zone-center. 
If the optical phonon branch has a maximum (minimum) value ujq at the zone-center, the 
cross section of the phonon-impurity scattering vanishes at u = ujq and increases for u"^ < uJq 
{uj"^ > uJq, respectively). Then the Raman line shape as a function of the frequency u shows 
asymmetry, being broader on the low (high) frequency side. The asymmetry is sensitive to 
the dimensions of the disorder |T^ . This important role of the " density-of-state" effect was 



noticed only recently for films subjected to laser heating |T^ and explain satisfactory the 
final line shape observed in this work. 

Our paper is organized as follows. To analyze the experimental data in detail, we dis- 
cuss the theory of probing the strain by Raman spectroscopy in Sec. II. In Sec. II-A, we 
recapitulate the equation of motion of optical phonons in the long-wave approximation and 
show that the effect of homogeneous strain gives rise to the splitting and shift of the optical 
phonon triplet. In Sec. II-B, we show how the static strain fluctuations involve the phonon 
scattering and result in increased phonon linewidths. Since the increase in linewidth can 
be larger than the natural phonon width in the anharmonic decay channels, we consider 
the effect of fluctuations self-consistently using Dyson's equation for the phonon Green's 
function. We obtain the inhomogeneous broadening and shift in terms of the strain corre- 
lation function. Supposing that the main strain effect is connected to the dislocations, a 
two-dimensional form of the correlation function is used. In Sec. II-C, we discuss the effect 
of surface roughness. In Sec. Ill, we summarize all experimental details and compare them 
with theory. In Subsec. III-A we discuss the results obtained from the infrared spectra 
collected as a function of SOL thickness for various SOI wafers. We show that they give 
evidence of residual permanent strain in the as-delivered wafers. Within the BOX, the strain 
relaxes more and more when the SOL thickness decreases. This is independent of the SOI 
wafer fabrication process as well as the thinning method. In Subsec. III-B, we show that 
this decrease of residual strain in the BOX correlates directly with an increase in the Raman 
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linewidth in the handle sihcon wafer. This is observed from micro-Raman spectra collected 
on the cleaved edges of wafers. Near the interface Raman spectra show that the threefold 
degeneracy of the A; ~ optical phonon is split indicating a large (local) strain. We deduce 
a noticeable inhomogeneous contribution to the linewidth, which extends down to several 
micrometers. Finally, we show that the surface SOL roughness is directly related to the 
strain fluctuations near the BOX/wafer interface. 

II. THEORY 

A. Summary of strain effect on optical phonons 

In Si, at the F point of the Brillouin zone the optical phonon displacements Ui{r,uj) 
(where i = x,y, z) belong to the threefold representation F25 of the diamond group 0\. The 
interaction of these phonons with the static strain eij{Y) corresponds to the third order terms 
(which are linear in ejj(r) and bilinear in ^^(r, 00) ) in the total energy expression. The first 
derivatives of the total energy with respect to ^^(r, 00) give then the equation of motion 

{{ujI + s^A - icuF(™*) - uj^)5,j + V.jiv)) Ujir, cu) = (1) 

where the zone-center optical phonon frequency ujq = 520 cm~^ for Si and the corresponding 
natural width F*^*"*^ is about 3 cm^^. The optical branches in cubic crystals are dispersive 
and, for instance, the dispersion of the x— polarized optical branch has the form: u;^ = 
^o — slk'^ — S2{ky + k'^). The neutron scattering experiments give an estimate of the dispersion 
parameters s^ ^ S2 ~ 0.8 ■ 10^ cm/s for Si. This shows that the dispersion term can be 
written in a more simple (isotropic) form s'^d'^/dxf in the vicinity of the branch maximum 
at the Brillouin- zone center. 

The interaction Vij{r) describes the changes in "spring constants" due to strain eij{r): 

''xx\'^) ■^xxxx^xxy'^ ) ' ■^xxyyl^yyy'^ ) ' ^zzy'^ )\ aUQ >'xy\'^) '^xyxy^ xyy'- ) ) \'^) 

where we take into account the fact that only three independent components \ijki exist in 
a cubic crystal. The other elements of l^jj(r) can be obtained by circular permutations of 



the indices in Eq. (^. It is evident that they are of order cJq. From experiments with 
homogeneous stress [O we find the coefficient values Xxxxx = P = —1-6 uJq, \xxvv = 1 = 



—2.2 Uq and Xxyxy = 2r = —1.3 Uq , where the generally accepted notation p, q, r is used. 

Because of the strain ffuctuations, all quantities have to be averaged over physically small 
distances (which are still large on a microscopic scale). The averaged strain {sij) conserves 
cubic symmetry in the plane parallel to the interface, but depends smoothly on the distance 
to the interface. In Fig. 1 this is chosen as the z— direction. We do not write this dependence 
in explicit form. The average components are such that: (Sxy) = {Sxz) = {^yz) = and 
{^xx) = {^yy)i which gives for the shifted phonon frequencies at the Brillouin-zone center 

iO^ik = 0) = (X'q + Xxxxx\^zz) + "^Xxxyyi^xx) (3) 

for the phonon polarization normal to the interface and 

'^x\^ = 0) = ^y{k = 0) = UJq + [Axxxx + Xxxyy)\£xx) + Xxxyyi^zz) (4) 

for the phonons polarized in the plane of the interface. This is nothing but the well-known 
phonon spliting into a 2;— singlet and an xy— doublet under the effect of a homogeneous 
strain |T^ - |T1 . 



B. Effect of strain fluctuations on the optical-phonon linewidth 

The contribution to the linewidth arises in the second-order perturbation theory for the 
strain components. It could be calculated using "the golden rule" of quantum mechanics. 
However, due to the singularity in the phonon density-of-states at u = ujq, this contribu- 
tion is strongly frequency-dependent and may be larger than the natural phonon linewidth. 
Therefore we calculate it self-consistently using a Green's function technique |18[. For the 



first-order Raman scattering, the corresponding cross section is bilinear in the phonon dis- 
placements. Then the Raman cross section 

^" oc ^^^^^^'^f"^,f":^:^' ImD^AKu^), (5) 



du^^'^dn^'') 1 - exp{-huj/kBT) 



is determined by the imaginary part of the phonon Green's function Daa'O^,^), where lo 
and k have the sense of frequency and momentum transfers, eL and e^'^^ are the polarization 
vectors of the incident and scattered photons, respectively, and ga/S'y is the /37— component of 
the Raman polarizability tensor for the a— phonon branch. It is easy to see that, for the r25 
representation, the tensor ga/3'y is completely antisymmetric and reduces to the component 
Qxyz- This means that, when the incident and scattered fields are polarized along the x and 
y axes, respectively, only 2;— polarized phonons can be excited. 

The phonon Green's function has to be averaged over the strain distribution. The calcu- 
lation has been performed in details in a previous work [Q and we shall simply summarize 
the results. The average phonon Green's function has the diagonal form: 

Djj(k, cu) = (j^5(k, uj) - s'^k^ - icuTjik, cu) - cj^)"^ , (6) 

where the unknown real functions Qj(k,uj) and rj(k, cj) obey the system of Dyson's equa- 
tions 

(7) 

The average interaction (Vjj) is defined in Eqs. (^) - @ and the correlation function in 
real space is 

W,Ur - r') = {6V,m{r)SV^,{r')) (8) 

with SVjmir) = Vjmir) - {Vjm)- 

Several points should be emphasized. First, one can see from Eqs. (|^) and (^) that 
the line center as a function of the frequency transfer u is determined by the equation 
r2|(k, u;) — s'^k'^ = uj^, where k is the momentum transfer. Here we can put k = 0, since the 
photon wave vector is small. Second, in the limiting case rj(k, uj) -^ 0, the bypass around the 
pole in the right-hand side of Eq. (^ gives the Dirac delta-function (5(fi^(q, u) — s^q^ — oj"^). 
We then arrive at the well-known golden rule and see that the strain fluctuations make a 



contribution into the linewidth only for uj < f2^, i. e. on the low- frequency side of the 
Raman line. The opposite situation occurs when the phonon branch has a minimum at the 
center of the Brillouin zone |]16l but, in both cases, the line shape shows asymmetry. 



C. Simplest model of two-dimensional disorder 

Equations (|^) can be simplified if we assume that the correlator Wjm{(i) has a constant 
value Wjmi^. = 0) in the region q < I/tq and vanishes elsewhere. On the other hand, the 
strain fluctuations may be considered as imperfections where Wjm{(\ = 0) is proportional 
to the concentration of imperfections and tq is the domain size in which a phonon interacts 
with the imperfection. The final form of Eqs. (^ depends on the dimensionality of the 



imperfections. In our previous work ||T4| we have concluded that the imperfections involved 
by interfaces are mainly two-dimensional, i. e. that they behave like dislocations. In that 
case, Tq is the dislocation core and the correlation function W{v) in Eq. (^ depends on the 
two-dimensional vector r^. Then Eqs. (|^) have the form: 

n|(^) -^,- {v,i) - 'iu{T,(^) - r<"">) 

where Bjm = Wjm{^± = 0)/47rs^ti;o5 and the sum can be taken over two split modes, which 
are the normally polarized singlet and the doublet parallel to the interface. Notice that the 
functions ^j{uj) and Tj^u) are independent of k (we omit k in the arguments) because of 



its smallness in comparison to the value y/uJoT/s ^ [tx /a)JT /ujq which is essential in the 
integral (|^). Again a is the lattice parameter. 

Let us estimate Bjm- According to Eqs. (|ID and (|]), V{r- = 0) = Ae ~ oJqS, where we 
omit the tensor indices. We have for the Fourier component (§) Wjmik = 0) ~ 2T:ujQ{5eYrQ 
and for the phonon-strain scattering probability 

B ~ ujl{r^6elsf (10) 
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where 6e is the strain fluctuation. 

If we assume that the strain fluctuations are mainly induced by dislocations, then 

B = cv^{k = 0)/4:7iuJos^. (11) 

In this expression, c is the dislocation concentrations, v{k = 0) = gr^UQ is the phonon- 
dislocations interaction where g is a dimensionless constant of the order of unity and Tq is 
the radius of a dislocation core. 

In the following we will solve numerically Eqs.(P) in ^j{uj) and Tj{uj) to flnd the values 
of Bjm and tq which give a flt of Eqs. (^ and @ to the experimental Raman spectra. 
Notice that the effect of the strain fluctuations on the line shape is qualitatively evident 
from Eqs.@. For instance, if there is only one phonon mode, its width T{uj) is governed by 
the equation 

,. , ,,,„R / ,,,2 _ 02|',,,^| j_ ^2/^2 ,,,2 _ o2^,,,^|\ 

(12) 



T(uj) - r(*"*) = -^ arctan ^-^- '-^ - arctan „, ) ' 



Let us consider the case of small tq (tq < s/vuJT). The terms in the right-hand side of 
Eq. (|T2|) compensate for each other on the high-frequency side [uJ^ > Q"^) of the Raman 
line, where the line has a nearly Lorentzian form with a width T ~ r'^*"*-'. For the low- 
frequency side (cu^ < Q"^), these two terms give a contribution to the linewidth which is 
proportional to the squared strain fluctuations. For the situation of a minimum in the 
phonon branch at the center of the Brillouin zone, the Raman line drops more rapidly on 
the low-frequency side. This effect has origin in the singularity of the phonon density- 
of-states at the branch extremum. In the opposite case, rg ^ s/\/uJT, the lineshape is 
symmetric but non-Lorentzian. 

Since T{uj) is tu— dependent, we must connect its value with the full width at half max- 
imum (FWHM) of the experimental spectra. The simplest way is to use T{uj) at the line 
peak, where the equation Q{uj) = cj is satisfied. In the case where F is independent of u, this 
definition gives the FWHM. The first term in the parenthesis Eq.([T^) remains only at the 
line peak. Finally, because the difference between the line positions Q{uj) and uq (including 
and ignoring the effect of disorder, respectively) is small, we obtain for F at the line peak 

11 



s^ 



Y - r(^"*) = B arctan ^ . (13) 

Notice that, in the hmiting case tq ^ s/a/cJT, Eq. ([l3|) coincides with the result of the 
perturbation theory T = r*^*"*^ + nB/2. In the opposite case arctan(s^/rQCt;or) = s'^/rQUJoT 
and we arrive at a quadratic equation which gives 

"n(mt) 



+ 



\ 






This shows that, in the general case, the effect of disorder which is represented by the 
parameter B, is not simply additive with the natural width r'^*"'*\ 

D. Surface roughness induced by dislocations 

Because of the interface strain and related dislocations, there must be a finite surface 
roughness contribution. Due to the very complicated nature of the interactions between the 
surface and dislocations it seems impossible to reach a high degree of sophistications in this 
problem. However, we can estimate the mean roughness using a qualitative approach. We 
suppose that the strain has large fluctuating components and calculate the mean square 
roughness. 

Any rough plane may be conceived as a set C,{s) = I]„&(s — s„) of randomly located 
"hills" at points s„, where s is a two-dimensional vector of the free surface plane. If we 
suppose that the hills originate from dislocations, the height b must be on the order of the 
Burgers vector Bq. We specify / as the size of the hills (i.e. the domain size where the 
function b{s — s„) has a nonzero value) and, by averaging over the distribution of points s„, 
we obtain the roughness correlation function 

C2(s - s') =< C(s)C(s') >= cY: I &(q) I' e'^'^^-^l 

where 6(q) ~ BqP is the Fourier component of the function 6(s) and c is the hill concentration 
at the surface. 
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The mean square roughness is given by ^2(8 = s') ~ cfeg^^ which is connected to the 
strain relaxation by the dislocation concentration c. Of course some initial roughness Co 
does exist, which is not strain-dependant. Taking both contributions into account, we write 
for the effect of strain on the surface roughness 

C^ = Co' + cbll'. 



This expression and Eqs. (]TT|), ( ]TBD enable us to compare the change in roughness with the 
phonon linewidth. We get 



r - r(-*) = ^(C^ - Co') arctan -^, (15) 



'1 
where A = g'^r^ujQ/ATTs'^blP. 

III. EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 

Most experimental details have already been given in |^ and will not be repeated here. 
To summarize briefly, three different (commercial) 4-inch wafers from SOITEC [^] were 
considered. One wafer was high dose SIM OX and two were Unibond. In all three cases, 
the nominal SOL thickness was 0.2 /im. In two cases (one Unibond wafer and the SIMOX 
wafer), a wet SO (Sacrificial Oxidation) process was used. This resulted in two series of 
samples denoted A and B, respectively. In order to evaluate the effect of the repetitive high 
temperature oxidation steps (at 1050 C) on the final results, a low temperature oxygen- 
assisted IBE (Ion Beam Etching) process was applied to the second Unibond wafer. This 
resulted in the third series of samples (samples C). 

A. Influence of the SOL thickness on the IR reflectivity spectra 



To control the SOL thickness |T^ infrared reflectivity spectra were systematically col- 
lected at room temperature with a Brucker IFS 66v Fourier Transform Spectrometer fitted 
with a DTGS detector. Measurements were made in the spectral range 400 to 7500 cm~^ 
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and a standard "near normal incidence" configuration (finite value of the angle of incidence 
close to 20°) was used. The direction of incidence was in the yz— plane of Fig. 1 and no 
attempt was made to polarize the incident light. 

Apart from the change in interference patterns which gives the final sample thicknesses, 
a strong effect was found in the 1000 - 1100 cm~^ range. This is the narrowing of the Si02- 
related feature displayed in Fig. 2 for samples A. It has nothing to do with the change in 
the interference pattern, but correlates with the change in SOL thickness. Notice that there 
is a large effect when reducing the SOL from 203 to 95 and 42 nm. Notice also that, starting 
from 42 nm, it seems to saturate. 

The 1070 - 1080 cm~^ doublet structure is well known. It manifests itself in every 
silicate-like material (not only in crystals, ||20|- |23|, but also in vitreous silica p4 



or, even, 



in sol-gel derived glasses [^) and corresponds with internal vibrations of the building Si04 
tetrahedra. It is very important to notice the stability of the corresponding frequency which 
does not change in going from j3— to a— quartz [^] or, even, when Si is replaced by Al or P 
in AIPO4, which is the ternary analog of quartz [^. In the literature it has been classified 
as an AS (Asymmetric Stretching) internal mode P3 , E^ in which one oxygen atom bonded 



to two neighboring Si moves parallel to the Si-Si direction. This oxygen motion is strongly 
infrared active and, despite the lack of long-range ordering, in many cases the associated 
LO frequency resolves around 1260 cm~^. Due to the finite value of the angle of incidence, 
this feature appears also in the present work. 

Concerning Fig. 2, a second point should be emphasized. In the specific case of a 
thin Si02 film on Si, one expects two equivalent in-plane x— and y— modes and only one 
out-of-plane z— mode. In other words, because we are dealing with a thin layer in which 
the threefold degeneracy is lifted, one expects to find a doublet structure with a 1 to 2 
intensity ratio. This is exactly what is found. Starting from one broad (unresolved) feature 
for nominal SOL thickness 200 nm, one goes rapidly to two distinct vibrational modes when 
the SOL thickness reduces below 40 nm. From the corresponding intensity ratio, we deduce 
that the 2;— mode appears lower. Since the average energy difference with respect to the 
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X, y— modes comes from the thermal strain, the strain is compressive. Making the length of 
the Si-0 bonds oriented in the xy— plane shorter, it shifts the x, y— modes to higher energy. 
Of course there are some finite deviations with respect to the average energies and, since we 
find evidence of a clear re-ordering (within the thin dielectric film) when thinning the SOL, 
we must conclude that the SOI system is not at all rigid but, instead, at the atomic scale 
constitutes a rather fiexible medium. 

Since the simplest way to shift the frequency of an internal mode is to change the 
bond length, our results suggest that we are dealing with oxide relaxation. To explain 
from the very beginning the narrowing of the Si-0 related feature, one must assume an 
initial (thermal strain induced) distortion which relaxes when the SOL thickness decreases. 
Starting from a broad feature (~100 cm^^) for the as-delivered material (broad distortion 
of the equilibrium bond length) the narrowing of the Si-0 bond signature demonstrates a 
reordering of the Si-0 distances as the SOL thickness decreases. This is true from the very 
beginning, when thinning the SOL from 200 to 100 nm, and does not depend on the type of 
material investigated. This observation has been checked on SIMOX (series B) and similar 
results were found. 

Finally, to check that this was not due to any parasitic effects of the high temperature 
SO steps, we repeated the thinning process using oxygen-assisted IBE (series C). Again, 
similar results were found. This lends support to the proposition that we are dealing with 
real and intrinsic properties: the experimental improvement found in the oxide relies on an 
internal stress relaxation mechanism, where the stress originates only from the difference in 
thermal expansion coefficients between Si and Si02 (~80% at room temperature). There is 
no artifact associated with the high temperature SO steps and, to the best of our knowledge, 
this provides the first clear experimental evidence that all SOI materials are plastic systems, 
soft enough to relax when the SOL thickness decreases. 

To estimate the relaxation we write the interatomic potential energy proportionally to 
r"", where n should be of the order of 5 for a typical oxide. Because the relative change of 
the vibration frequency is 6uj/uj = {n + 2)6r/2r and the linewidth of the refiectivity signal 
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changes by around 50 cm ^ when thinning the SOL, we deduce that the average change in 
Si-0 bond length is a reduction of about 1%. 

B. Raman investigation of nominal and thinned SOI 

We have mainly shown that, because of internal stress relaxation, the BOX properties 
improve when thinning the SOL. Now we consider what happens in the silicon overlayer 
and wafer. First we investigated in details the Si handle-wafer by using micro-Raman spec- 
troscopy. We use a Jobin-Yvon ISA T 64000 spectrometer fitted with an Olympus micro- 
scope and a cooled CCD detector. The 488 nm line of a Spectra-Physics argon-ion laser was 
used as the excitation frequency. To achieve in-depth resolution, a transverse backscattering 
configuration was used [|1^] with both the incident and scattered light propagating along the 



x— direction of Fig.l. Because of the large experimental aperture we probe, both, y and z 
polarized phonons. 

Typical results are shown in Figs. 3 to 6 for the series of samples A (Unibond). Fig. 3 
corresponds to nominal material with SOL thickness 203 nm (standard value). Fig. 4 with 
42 nm (thinned). Fig. 5 with 11 nm (thinned) and, finally. Fig. 6 with 5 nm. Given a specific 
SOL thickness, the different spectra correspond with the different laser spot distance from 
the wafer surface. To quantify in full details the different contributions, they have been 
fitted using the theory of phonon interaction with static strain fiuctuations presented in 
Sec. II. The results are shown as solid lines on the different figures. The complete series of 
parameters is given in Table 1. 

The fits with Eq. (^ and Eq. (^ were made simultaneously for all spectra. The best 
value of the parameter tquoq/ s = 16 gives a reasonable value for the correlation radius (or the 
dislocation core) tq = 1.3A. Eqs. (^ were solved numerically to obtain the matrix elements 
Bij, which determine the inhomogeneous linewidth and the shift of the lines. One example 
of the solution (functions Tj{uj) and ^j{uj) — uq— {Vjj)/2) is presented in the right panel of 
Fig. 4 for the z— phonon {j = z). The resulting values (i.e. the matrix elements B^^ and 
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Bzz) are given in Table I for samples A. 

Every time, focusing far from the interface (a typically value is around 5 fim. for the 
complete series of samples) one finds an identical (standard) "bulk" Raman linewidth for 
the silicon handle wafers. In the following this is taken as "reference material" and the line 
position for the reference material (520 cm~^) is shown by the vertical dash-dotted line as 
a guide for the eye. Moving toward the surface (first to 2 /xm and, then, to 1 /im) has two 
consequences: i) one shifts significantly the phonon frequency to lower energy and ii) one 
broadens considerably the width of the Raman peak, especially in the cases of the thinner 
SOL. This is evident starting from 42 nm. 

A first contribution to the phonon lineshift towards the low frequencies comes as a result 
of the expansion of the Si lattice near the Si/Si02 interface ||12|, H^. This homogeneous shift 



is described by Eqs. (0) - (|) but, because of the small energy difference, the strain-induced 
splitting into a singlet and a doublet states resolves only for the thinnest SOL thicknesses 
(5 and 11 nm, respectively: see Table 1). 

The second contribution comes from the strain relaxation mechanism which manifests it- 
self experimentally when probing at different positions along the z— direction. It is frequency 
dependent and gives the inhomogeneous linewidth and shift contributions displayed in the 
right panel of Fig. 4. Notice that, when focussing on the active part of the SOI wafers 
(upper part of the handle Si wafer) the inhomogeneous linewidth and shift contributions 
at the line center become very large. They range from 1.1 cm~^ (width) and +0.14 cm~^ 
(shift) for the 42 nm SOL sample to 6.1 cm~^ and +0.36 cm~^, respectively, for the 5 nm 
SOL sample (see again Table I). The inhomogeneous contribution to the line-shift remains 
however smaller than the homogeneous ones: -0.9 cm~^ for the 42 nm SOL sample and -1.2 
cm~^ for the 5 nm SOL sample. 

Similar results are found for SIMOX. They are displayed in Fig. 7 for the SOL thickness 
50 nm. The natural width is slightly larger (3.4 instead of 3.1 cm~^ for Unibox) but the 
correlation radius has the same value tqujo/ s = 16. The main distinction is that the depth of 
strain relaxation is larger: there is still a noticeable difference in the linewidths at distances 
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5 and 20 nm to the interface (see Table II). Finally, some homogeneity at the interface is 
found. Two spectra (last rows of Table II) collected from different points at the interface 
show practically the same characteristics. However other regions were detected, where more 
dramatic differences do exist. We shall come back to this point in Sec. IV (Conclusions). 

To summarize this section, we have found from our analysis of the linewidth variation 
versus distance that all strain-induced effects (like the homogeneous and inhomogeneous 
shifts and broadening) have a maximum near the Si/Si02 interfaces. Because they depend 
directly on the SOL thickness, they should affect primarily the thinned SOLs. The surprising 
effect is that they affect also the upper part of the handle wafers and relax through a depth 
on the order of several /im. 

C. Surface roughness and Raman ^vidth 

From the Raman investigations, we have found that the morphology of the upper part of 
every (nominal and thinned) SOI wafer experiences a complex stress distribution. Moreover, 
the strain in the wafer near the BOX/wafer interface becomes larger when thinning the SOL. 
In this section we show that independent evidence can be found from the consideration of 
root mean square roughness of the SOL surface. We reported already that there is a strong 
degradation of the surface when the SOL thickness decreases p. Typically, decreasing 
the SOL thickness from 203 to 5 nm, one increases the mean square AFM (Atomic Force 
Microscopy) roughness from 1.5 to 4 A RMS. 

To show that the Raman width and mean surface roughness have the same origin (inho- 
mogeneous strain and dislocations) we show in Fig. 8 a correlation of the Raman linewidths 
(cm~^) obtained at the surface of wafers A with different SOL thickness with AFM data 
(A). The solid lines are plots of Eq. (|15|). The experimental values are shown by dots. A 
single parameter A = 2.1 cm~^A~^ was chosen to fit the values of the 2;— phonon width and 
roughness in the case of the SOL thickness of 5 nm. Notice that the doublet splitting resolves 
only for the 5 and 11 nm SOL thickness. The value A = 0.28 cm~^A~^ for the x— phonon 



was obtained in the comparison of the interaction constants Bxx and Bzz from Table I. From 
these data, it is interesting to estimate the parameter / of the surface roughness using this 
value of A (ro is obtained in this work and s is known). Taking into account the expression 
of A given after (^) and putting g = I and 60 = 1^, we get a very reasonable value / = 8SA 
for the average diameter of the hills. The main difference from one thinning process to the 
next one is in the initial roughness. With this respect, all results obtained using SO (1.2 A, 
series A) appear better than the one obtained with IBE (2 A, series C). 

IV. CONCLUSIONS 

It has been known for a very long time that any direct bonding of two oxidized wafers (or 
any oxygen-ion implantation followed by a high temperature annealing step, for instance) 
introduces a finite stress at the Si/Si02 interface. This comes directly from the difference 
in thermal expansion coefficients between Si and Si02 (2.6 ■ 10^^ K~^ and 0.56 ■ 10^^ K~^, 
respectively 0). Because this constitutes a hard point of the technology, low temperature 
bonding using low viscosity oxides like BPSG (BoroPhosphoroSilicate Glass) or SOG (Spin 
on Glass) are under active development. In some cases they will offer alternative solutions. 
The problem of process- induced defects is totally different. To the best of our knowledge, it 
is much less documented. For instance there is no clear report of a critical SOL thickness, 
similar to the critical layer thickness encountered in hetero-epitaxy. Above such a critical 
thickness, the SOL material would be homogeneously strained but stable. It would strongly 
relax around that thickness, with a maximum in inhomogeneous behavior. Finally, it would 
be heavily strained, but more homogeneous and again stable below that value. Our experi- 
mental results suggest that such a critical thickness does exist. It is about 40 to 50 nm and 
corresponds with the experimental situations illustrated in Fig. 2 (BOX), Fig. 9 (Unibond) 
or Fig. 10 (SIMOX). In all three cases we find a critical regime where, simultaneously, the 
BOX achieves a better relaxation while there is a very large inhomogeneity in the SOL. 
The spectra displayed in Figs. 9 and 10 (solid lines are again theoretical curves) could be 
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understood in the terms of phonon splitting of the order of 7 to 15 cm~^ due to a very high 
local strain in some finite parts (or domains) at the interface. 

To summarize, we have evidenced in this work that SOI is not a perfectly strain-relaxed 
system. It behaves more like a balanced-strain structure schematically drawn in Fig. 11. 
The interesting point is that: i) compressive strain (about 1 % in absolute value) is always 
present in the buried oxide of as-delivered SOI wafers. This strain diminishes when the SOL 
thickness decreases; ii) on contrary, tensile strain exists in the silicon overlayer. Dislocations 
near the interface release the strains and manifest themselves in final surface roughness (up 
to 4 A RMS). The dislocations can be considered as the origin of the strain fluctuations 
which result in optical phonon scattering; finally, iii) tensile strain and strain fluctuations 
(again on the order of 10^^) also exist in the bulk of the underlying wafer. Every time, the 
strain originates from the Si/Si02 interfaces and comes from the difference in bond length 
and thermal expansion coefficient between the two materials. Since both the Si-overlayer 
and the bulk Si wafer keep the buried oxide stressed, when thinning the SOL the strain 
decreases in the buried oxide and increases in the handle wafer. The overall strain relaxes 
very slowly and extends from the BOX/wafer interface over a few micrometers range. This 
is in contradiction to the most common belief that the SOL/BOX/underlying handle wafers 
are made of very stable and strain-free materials. As a matter of fact, SOI wafers experience 
more equilibrium conditions which depend on every stage of the technology road. 
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Figure captions: 

Fig. 1. Schematic drawing of tlie Raman back-scattering geometry used in this work to 
investigate the strain relaxation in the silicon-on-insulator (SOI) system. The strain relax- 
ation is probed by displacing the laser spot in the ^-direction. 

Fig. 2. Infrared reflectivity spectra collected on SOI (Unibond) material for different SOL 
thickness (see labels on the curves). 

Figs. 3. Experimental Raman spectra collected on Si-wafers with SOL thickness 203 nm. 
The SOL is separated by a 400 nm oxide film from the handle wafer. The distance from the 
interface (5, 2, 1, and /im ) of the laser spot is indicated on the curves. The solid lines are 
fits to Eqs. (§P - (|3)- The results of the fits are given in Table 1. The natural width is 3.1 
cm^^ for all lines, the correlation radius tqUq/ s = 16. 

Fig. 4. Same as Fig. 3 for 42 nm Si-overlayer thickness. In right panel, the inhomogeneous 
shift and linewidth are shown as a function of the frequency. The position of the line peak in 
fully relaxed material is shown as dash-dotted line. The inhomogeneous shift (right panel) 
is positive, whereas the total shift (left panel) resulting from the strain is negative. 
Fig. 5. Same as Fig. 3 for 11 nm SOL thickness. 
Fig. 6. Same as Fig. 3 for 5 nm SOL thickness. 

Fig. 7. Same as Figs. 3-5 but now for a SIMOX wafer with 50 nm Si-overlayer thickness. 
Results of the fits are given in Table II. The natural width is 3.4 cm~^ for all lines, the 
correlation radius tqUq/s = 16. 

Fig. 8. Linewidths of the two components of the optical phonon triplet in Si near the Si/Si02 
interface as a function of the root mean square (RMS) roughness of the SOL surface. The 
lines are the result of the theoretical model (see Eq. (|15|) ). The points give the width found 
from the fit of spectra to the theory. The corresponding RMS value is extracted from AFM 
measurements for the various thickness (203, 95, 42, 11, and 5 nm) of the Si-overlayer. 
Fig. 9. Experimental Raman spectra (for Unibond wafers A) collected for different points 
at the BOX/wafer interface (upper line: normal situation; bottom lines: specific points dis- 
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cussed in the text). The thickness of the Si-overlayer is 42 nm. 

Fig. 10. Full series of Raman spectra (for SIMOX) collected up to 5 ;um from a specific 

point at the BOX/wafers interface. The thickness of the Si-overlayer is 50 nm. 

Fig. 11. Schematic drawing of the in-plane strain in the Si-overlayer, the buried oxide and 

the handle silicon wafer for two different SOL thickness: nominal value 200 nm (solid line) 

and thinned SOL (dashed line). 
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TABLES 
TABLE L Zone-center optical phonon frequencies and corresponding linewidths (all values are 

in cm"^) in the Unibond samples (with various top-layer thickness) in relation to the distance z 

from the interface. The results refer to the fits with the present theory. The adjusted coupling 

parameters (the last column) B^x^^zz are probabilities for the phonon intrabranch scattering by 

strain fluctuations. The splitting of triplet r25 into the doublet iyix) and the singlet (w^) is resolved 

for samples with 5 nm and 11 nm top-layer thicknesses. 



26 



Sample 


Distance from 


Lines Fgs 


Linewidths 


Interactions 




interface, /im 


w^; LOx 


r.; r. 


Bzz ; Bxx 


203 nm 


5.0 


520.0 


3.10 


0. 




2.0 


519.76 


3.21 


0.20 




1.0 


519.28 


3.32 


0.39 







519.18 


3.61 


0.97 


95 nm 


5.0 


520.0 


3.10 


0. 




2.0 


520.0 


3.10 


0. 




1.0 


520.0 


3.32 


0.38 




0. 


519.28 


3.69 


1.16 


42 nm 


5.0 


520.0 


3.10 


0. 




2.0 


520.0 


3.21 


0.19 




1.0 


519.18; 516.13 


3.69; 3.26 


1.16; 0.29 




0. 


519.17; 519.08 


4.25; 3.42 


2.51; 0.6 


11 nm 


5.0 


520.0 


3.10 


0. 




2.0 


520.18; 520.15 


3.32; 3.1 


0.38; 0. 




1.0 


520.27; 520.22 


3.86; 3.42 


1.54; 0.58 




0. 


519.55; 519.13 


6.23; 3.70 


9.6; 1.2 


5 nm 


5.0 


520.0 


3.10 


0. 




2.0 


519.73; 519.70 


3.32; 3.1 


0.39; 0. 




1.0 


519.26; 519.27 


6.44; 3.61 


10.6; 0.96 




0. 


519.14; 519.00 


9.44; 4.65 


29.0; 3.67 
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TABLE II. Same as Table I but for a SIMOX sample with 50 nm top-layer thickness. The 
splitting of triplet was not resolved 



Distance from 


Lines r25 


Widths 


Interactions 


interface, fim 


^x 


r. 


Bxx 


20.0 


520.0 


3.40 





5.0 


520.01 


3.50 


0.20 


2.0 


520.02 


3.60 


0.38 


1. 


519.76 


3.70 


0.58 


0. 


519.50 


3.79 


0.77 


0. 


519.30 


3.79 


0.77 
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